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Aims. We present our analysis of the multiwavelength photometric & spectroscopic observations of GRB 060210 and discuss the 
results in the overall context of current GRB models. 

Methods. All available optical data underwent a simultaneous temporal fit, while X-ray and y-ray observations were analysed tempo- 
rally & spectrally. The results were compared to each other and to possible GRB models. 

Results. The X-ray afterglow is best described by a smoothly broken power-law with a break at 7.4 hours. The late optical afterglow 
has a well constrained single power-law index which has a value between the two X-ray indices, though it does agree with a single 
power-law fit to the X-ray. An evolution of the hardness of the high-energy emission is demonstrated and we imply a minimum host 
extinction from a comparison of the extrapolated X-ray flux to that measured in the optical. 

Conclusions. We find that the flaring y-ray and X-ray emission is likely due to internal shocks while the flat optical light curve at 
that time is due to the external shock. The late afterglow is best explained by a cooling break between the optical and X-rays and 
continued central engine activity up to the time of the break. The required collimation corrected energy of ~ 2 x 10 52 erg, while at the 
high end of the known energy distribution, is not unprecedented. 

Key words. Gamma rays: bursts - X-rays: individuals: GRB 060210 - Dust, extinction - Radiation mechanisms: non-thermal 



1. Introduction 



Since the launch of the Swift satellite (Gehrels et al. 2004b we 
have been allowed observations of Gamma Ray Bursts (GRBs) 
at very early times, occasionally observing in X-rays and optical 
light even as the prompt y-ray emission is ongoing. Gradually, 
evidence has emerged for a fairly canonical behaviour of y-ray 
and X-ray light curves in which the X-rays decay fairly rapidly 
after the prompt emission and then decay more slowly for a 
while before becoming steeper again dNousek et al. 2 006). This 
behaviour is seen in many of the bursts, but a significant minority 
decay more gradually from early times. These differences may 
be due to the different relative strengths of emission components 
due to internal and external processes dO'Brien et al. 2 006). 
Similar complexity is seen in the optical, where some bursts 
show optical behaviour which does not exactly mimic the 
X-ray, again suggesting a contribution from several emission 
components. Here we present multiwavelength observations of 
GRB 060210, covering the early to late-time emission and com- 
pare the data to the expected canonical behaviour and the blast 
wave afterglow model (e.g. Zhang & Meszaros 2004). 

GRB 060210 is a long burst, d etected by the Burst 
Alert Telescope (BAT) on board Swift ( |Barthelmy et al. 2005| ) 
on February 10 th, 2006 at 04:5 8:50 UT. The X-Ray 
Telescope (XRT; IBurrows et al. 20031 ). which started observ- 



ing the region 95 seconds after the BAT trigger, identi- 
fied the X-ray counterpart position to within an error of 
5.4" dBeardmore et al. 20061). While t he Ultraviolet/Optical 
Telescope (UVOT; |Roming et al. 2005) was unable to iden- 
tify an optical counterpart, the robotic Palomar 60" telescope 
observed the R-band counterpart at approximately 5.5 min 
post trigger dFox & Cen ko 2006). The redshift was estimated 
as z = 3.91 ( |Cucchiara, Fox, & Berger 2006]) and a host 
galaxy has been proposed based on near infrared observations 
dHearty et al. 2006) . No radio emission has been observed down 
to a 3cr limit of 72 /Jy at 8.64 GHz dFrail 20061) . 



Send offprint requests to: P.A. Curran (pcurran@science.uva.nl) 



In this paper we present our multi-wavelength analysis of 
the burst from y-ray, through X-ray to optical wavelengths. 
Throughout, we use the convention that a power-law flux is 
given as F v oc r a v~P where a is the temporal decay index and 
/3 is the spectral index. The power-law of spectra is given as 
dN/dE oc E~ r where F is the photon index and is related to spec- 
tral index as T = + 1. In ^2] we introduce our observations and 
reduction methods of each wavelength regime. In |3]we present 
the results of our spectral and temporal analysis, while in ^4] we 
discuss these results in the overall context of the burst. 
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Fig.l. Full BAT light curve (15 - 350 keV) of GRB 060210, 
clearly showing the main peak and post-trigger activity, as well 
as the pre-trigger activity at ~ To - 230 s. 

2. Observations 

2. 1 . y-ray 

The Swift-BKT triggered and localised GRB 060210 at T = 
04:58:50 UT and promptly distributed the coordinates via the 
Gamma-ray burst Coordinates Network (GCN). Data for the 
burst were obtained in event mode, covering To - 300 s to 
To + 300 s, giving spectral coverage from 15 - 350 keV and time 
resolution of 64 ms. 

BAT event data were reduced and analysed within FTOOLS 
v6.0.4. Light Curves and spectra were extracted by the standard 
methods, and standard corrections applied (O'Br ien et al. 2 006). 
From the light curve (figure Q} the main period of activity is 
identified as being from To - 75 s to To + 20 s with the bright- 
est peak at To- However significant emission was recorded from 
To - 230 s to the last peak at To + 200 s, which represents a sig- 
nificant fraction of the total recorded event data. The 7^0 in the 
BAT range ( 15-350 keV) is calculated as 220 + 70s, covering the 
period from ~ To — 149 s to ~ To + 70 s. 

2.2. X-ray 

The XRT started observing at 05:00:25 UT, 95 seconds after the 
initial BAT trigger, and found an uncatalogued, variable X-ray 
source which was identified as the afterglow. XRT observed up 
to March 3, or approximately 2 x 10 seconds post-burst, by 
which time only an upper limit detection was obtainable from the 
data. Data were initially collected in Windowed Timing (WT) 
mode (IHill et al. 20041 IHill et al. 20 05). which gives no posi- 
tional information but high temporal resolution. After the initial 
~ 500 s of WT mode, data were collected in Photon Counting 
(PC) mode. 

The XRT data were initially processed with the FTOOL, 
xrtpipeline (vO.9.9) with wtbiasdiff applied to correct for 
potential WT bias-row subtraction problems, caused by bright 
Earth or CCD temperature variations. Background subtracted 
0.3-10keV light curves were extracted from the cleaned event 
lists with a minimum of 20 counts/bin. Likewise, source and 
background spectra were extracted for analysis with Xspec. The 
first orbit of PC mode data were corrected for pile-up as detailed 
in |Vaughan et al. (2 006 ). A correction was also made, using 



Table 1. FTN optical observations of GRB 060210 calibrated to 
the Cousins photometric system, uncorrected for Galactic ex- 
tinction of E^B-v) - 0.093. Magnitude errors are at lcr level. 
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65.58 


31.33 
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22.10 fl 


5.85 


0.1667 
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18.50 ± 0.20 


6.20 


0.1667 
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18.36 ± 0.15 


6.55 


0.1667 
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18.36 ±0.15 


11.85 


0.5 
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18.51 ± 0.12 


16.35 


1.0 
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18.87 ± 0.10 


23.20 


2.0 


R 


19.43 ± 0.10 


33.04 


3.0 
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19.74 ± 0.10 


43.25 


2.0 


R 


20.14 ± 0.15 


53.18 


3.0 


R 


19.98 ±0.15 


72.00 


2.0 


R 


20.57 ± 0.20 


88.84 


5.0 


R 


20.8 ± 0.4 


115.26 


9.0 


R 


21.2 ±0.3 


167.0 


30.0 
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22.4 ± 0.5 


9.43 


0.1667 
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16.94 ± 0.08 


13.19 


0.5 
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17.20 ± 0.05 


18.11 


1.0 
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17.66 ± 0.05 


25.95 


2.0 


I 


18.25 ± 0.06 


36.98 


3.0 




18.56 ± 0.06 


46.00 


2.0 




18.80 ± 0.08 


57.10 


3.0 




19.18 ±0.08 


70.12 


0.1667 




19.35 ± 0.25 


73.73 


0.5 




19.43 ± 0.20 


78.73 


1.0 




19.55 ± 0.20 


86.41 


2.0 




19.85 ± 0.30 


97.26 


3.0 




20.0 ± 0.3 


106.13 


2.0 




20.2 ± 0.3 


116.97 


3.0 




20.3 ± 0.3 



a 3cr upper limit 



xrtexpomap, for fractional exposure loss due to bad columns 
on the CCD, which arose after damage caused by a microme- 
teroid strike ( |Abbey et al. 2006[ ). 

Spectral fits are used to convert count rates of both WT & PC 
mode data to fluxes, and to extrapolate the BAT light curve to the 
XRT energy range dO'Brien et al. 20 06 , figure|2]l. Initial analysis 
of the early time part of this light curve shows two strong flares 
at ~200s and ~380s, after an initial decay that seems to be a 
continuation of a flare at ~100 s as observed in y-rays. 

2.3. Optical 

Optical observations in B, R and i' bands were obtained by the 
robotic, 2.0 m Faulkes Telescope North (FTN) at Haleakala on 
Maui, Hawaii; FTN is the sister instrument of The Liverpool 
Telescope dGuidorzi et al. 20061) . At 05:04:36 UT, approxi- 
mately 5.8 min post-burst, the telescope triggered its automatic 
1 hour-sequence and to ensure good time coverage of the light 
curve, a second 1 hour sequence was manually triggered at 
06:05:27 UT, 66 min post-burst. The previously reported opti- 
cal transient (Fox & Cenko 2006 ) was clearly visible in R and i' 
but not in B (Table [l}. This is not unexpected, as the high red- 
shift of the burst and the associated extinction due the Lyman 
forest causes significant dimming in the B-band while i' remains 
unaffected. 



P.A. Curran et al.: The prompt to late-time multiwavelength analysis of GRB 060210 



3 



> 



I 

CO 

d 



3 




100 1000 10* 
time since trigger (s) 

Fig. 2. High-energy light curve of GRB 060210 composed of the 
XRT data (red crosses) and the BAT data (blue boxes) extrapo- 
lated to the XRT energy range. The dashed line shows a single 
power-law fit to the late afterglow with a = 1.1 1 ± 0.02, but the 
solid line showing a smoothly broken power-law with indices 
ai = 0.89 + 0.05 & a 2 = 1.35 ± 0.06 is a significantly better fit. 



The data were analysed using the IRAF packagefl, 
wherein reduction & differential photometry were carried 
out. Landolt photometric standards taken on the night were 
used to calibrate the magnitude of the optical counterpart. 
The i' magnitudes were converted to the corresponding I- 
band magnitudes using the observed filter transformations of 
Smith et al. (2002), which for normal colours are in agree- 
ment with those of Fukugita et al. (19%| . The magnitudes 
are uncorrected for Galactic extinction of E (B - V ) = 0.093 
( |Schle gel, Finkbeiner, & Davis 1998| l. Astrometric calibration 
was carried out with respect to the NOMAD catalogue^] to re- 
fine the position of the burst to within 0.2" to be 03:50:57.36, 
+27:01:34.4 (J2000). 



3. Results 

We have split our discussion of the emission into two phases, 
covering what we refer to as the early emission and the late after- 
glow. We take the early emission as being the y-ray and the flar- 
ing X-ray emission at times < 600 seconds, which follows from 
the X-ray light curve (figure [2]). We also include up to ~ 600 
seconds of optical observations, those not displaying a power- 
law decay, as being early emission. The late afterglow, we take 
as the smooth X-ray and optical light curves, obeying power-law 
decay. All uncertainties of light curve analysis, spectral fits and 
other are quoted at the 90% confidence level. 

3.1. Early high-energy emission 

The prompt emission as measured by BAT shows a main peak 
which decays smoothly until ~10s (figure From ~10s to 
~600s, GRB 060210 is observed to display flaring in both the 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

2 http ://w w w. nof s . navy, mil/nomad/ 



y-rays and X-rays. One flare, at ~100s, is observed by both 
BAT & XRT. This flaring behaviour seems to be superimposed 
on an underlying power-law decay due to the afterglow. The 
time-averaged BAT spectrum (15 - 150keV) over Tgo was fit 
with an unabsorbed power-law in Xspec and a photon index, 
r = 1.55 ± 0.09 was found (xl = 0.84, 56 d.o.f.). This spectral 
fit corresponds to a fluence of (6.0+q') x 10~ 6 erg cirT 2 . 

When the XRT WT-mode spectrum, from 103s - 614s 
is fit with an absorbed power-law, we find a photon index 
of T = 2.09 + 0.04 and a column density of N H = (16 + 
1) x 10 20 cirT 2 (xl = 1.020, 316 d.o.f.). This is significantly 
higher than the Galactic value of A^ H = 8.5 x 10 20 cirT 2 
( |Dickey & Lockman 1990] >. This excess extinction may be ex- 
plained either by rest frame extinction or by a broken power-law. 
Fitting an absorbed broken power-law gives a column density 
consistent with that of the Galactic value, so we fix the param- 
eter to be the Galactic value. This leads to a photon index of 
r = 1.96 ± 0.03 above a break energy, E hteak = 0.71 + 0.07 keV 
and an index of T = 0.6 + 0.3 Of 2 = 0.923, 315 d.o.f.) be- 
low the break. Alternatively, assuming solar abundances, a rest 



frame column density of A^h = (2.3 + 0.3) x 10 cm in com- 
bination with the Galactic value gives an equally good fit (xl - 
0.936, 316 d.o.f.) and implies a photo n index, T = 2.03 + 0.03. 
Assuming LMC and SMC abundances (Pei 1992) in the host im- 
plies N H = (5.7 + 0.7) x 10 22 cirT 2 (xl = 0.948, 316 d.o.f.) 
and Nh = (10.3 + 1.2) x 10 22 cirT 2 (xl = 0.950, 316 d.o.f.) 
respectively, while leaving the photon index unchanged within 
uncertainties. Though it is not possible to favour one model over 
the other, we see that the photon index, F ~ 2.0 is marginally 
different to that of the earlier y-ray spectrum. A hardness ratio 
(l-10keV/0.3-l keV) plot of the X-ray emission (figure[3]l con- 
firms spectral evolution, though the y-rays are too faint to obtain 
a similar plot. This evolution is similar to that observed in the 
y-rays for a number of BATSE bursts dBhat et al7 l994). 



Table 2. Rest frame column densities as found in combination 
with Galactic extinction. 



(xlO 22 cm" 2 ) 


SMC 


10.3 ± 1.2 


LMC 


5.7 ± 0.7 


Galactic 


2.3 ± 0.3 



Using the average photon index of the BAT & XRT spectra, 
we get a conversion from BAT count rate to unabsorbed flux 
in the 0.3 - lOkeV range of 4.38 x 10~ 7 ergcrrT 2 cts -1 . Using 
the same fit of the XRT spectrum, yields a conversion of 4.83 x 
10~ n erg cm -2 cts -1 which we apply to the WT mode data. 



3.2. Late X-ray afterglow 

The late X-ray afterglow was observed from ~lhour to ~23 
days, at which stage only an upper limit determination was pos- 
sible. The light curve for this period is quite smooth, with no 
obvious signs of flares or bumps, and a hardness ratio plot of 
the afterglow shows no signs of evolution. The spectrum gives 
results comparable with the prompt spectrum ( 33.11 ). including 
an excess extinction over that of the Galactic value. A broken 
power-law with Galactic absorption suggests a photon index of 
F = 2. 1 3 + 0.06 above a break at £ break = 1 .04 ± 0. 1 5 ke V and an 
index of F = 1.4 + 0.2 (xl = 1.021, 225 d.o.f.) below the break. 
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Fig. 3. (top) Light Curve of the Windowed Timing mode data 
of GRB 060210 and (bottom) the hardness ratio (l-10keV/0.3- 
1 keV) plot which clearly displays an evolution coupled to the 
count rate. 
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Fig. 4. Faulkes Telescope North R and I band light curves, sup- 
plemented by previously published data, shifted to the FTN R- 
band offset. Fit shown is a 3-index smoothly broken power-law. 



Assuming a single power-law with Galactic absorption, the rest 
frame absorption assuming solar abundances is consistent with 
that found during the prompt emission and hence is fixed at that 
value. The photon index is then T = 2.14 + 0.03 (xl = 1.015,227 
d.o.f.). It is not possible to favour one model over the other but 
we can assume a corresponding spectral index,/? = 1.14 + 0.03, 
at least at energies > l.OkeV. Using the single power-law fit 
with rest frame absorption yields a PC mode, count rate to flux 
conversion of 5 . 1 2 x 1 0~ 1 1 erg cm -2 cts" 1 . 

Fitting a power-law to the temporal decay we find a = 
1.11 ± 0.02 (xl = 1.69, 144 d.o.f.), though an F-Test shows 
that a smoothly broken power-law with a\ - 0.89 + 0.05, 
a 2 = 1.35 + 0.06 and a break at f break = 7.4+ 2 ;' hr (xl = 0-954, 
142 d.o.f.) is a significantly better fit, in agreement with the anal- 
ysis of |Dai & Stanek (2006} . Extrapolating the power-law fit to 
early times, we see that it does not match the early emission, 
which implies a break at ~ 10 3 s. Doing likewise for the broken 
power-law fit, we see that it matches quite well with the troughs 
of the flares and the extrapolated BAT emission (figure[2]i. To in- 
vestigate any possible spectral difference, we split the spectrum 
into pre- and post-break sections but the resultant fits were con- 
sistent with each other and the combined fit. This implies that 
any difference in spectral indices must be smaller than the er- 
rors, 60 < 0.05. 

3.3. Optical emission 

Fitting our optical light curves of the late afterglow (T > 500 s) 
with a single power-law, we find ctr = 1 .07 ± 0. 1 1 (xl = 1 -07) 
and ai = 1.19 ± 0.05 (xl = 0.86). To make full use of the avail- 
able data we simultaneously fit our data and previously publised 
data, at times > 540 s, over 5 bands: FTN R, FTN I, M DM R, 
KAIT I, & KAIT unfiltered dLi 20061 IStahek et al. 20071) . To do 
this we assume that each has the same temporal decay but we 
make no assumptions regarding the relative offsets - thus elim- 
inating any error in absolute magnitude estimates. This fit uses 
the simulated annealing method (§ 10.9 of lPress et al. 19921 and 
references therein) to minimise the combined^ 2 of the data sets 
and hence find the best fit parameters, while a Monte Carlo anal- 
ysis with synthetic data sets is carried out to estimate the errors. 



From this we find an optical temporal decay, a opt = 1.15 + 0.04 
(xl = 1.12, 38 d.o.f.) and magnitude offsets, which are the model 
magnitudes evaluated at 1 second, as detailed in table[3] This de- 
cay is shallower than the a ~ 1.3 found by Stanek et al. (2007} . 

Table 3. Simultaneous temporal fit of optical data, where Mag is 
the model magnitude evaluated at 1 second. The common tem- 
poral decay index, a opt is fit as 1.15 ± 0.04. 

Band Mag 

FTN R 10.22 ± 0.34 

FTN I 8.92 ± 0.30 

KAIT I 9.28 ± 0.30 

KAIT unfiltered 9.81 ± 0.30 

MDM R 10.11 ±0.34 



Shifting all the data to the FTN R magnitude offset confirms 
that there is no structure above noise in the optical light curve 
(figure|4|i. It is possible that the light curve is breaking at ~ 5000 
seconds, but there are not enough data to support this. Fitting 
a smoothly broken power-law with three indices to this shifted 
data, we find that the two break times are at 3 10+8 s and 540+6 s 
and the indices are a\ = 0.10 + 0.10, a 2 = -0.90 + 0.35, ai = 
1.175 ± 0.016 (xl = 0.89, 47 d.o.f.). 



4. Discussion 

It has been claimed ( Stanek et al. 2007) that this burst closely 
resembles GRB 050801 (Ry koffet al. 2006) > which is true inso- 
far as both display an initially flat optical light curve. However 
GRB 060210 displays a peak before decaying as a power law, 
while GRB 050801 shows no such peak. Furthermore, the X-ray 
behaviour of GRB 050801 seems to be consistent with the opti- 
cal, but in GRB 060210 they seem to be unconnected. This sug- 
gests that in GRB 050801 the X-ray and optical emission is orig- 
inating in the same region, while this is not true for GRB 060210. 
In the case of GRB 060210 the early flaring X-ray & y-ray emis- 
sion implies internal shocks, indicative of prolonged engine ac- 



P.A. Curran et al.: The prompt to late-time multiwavelength analysis of GRB 060210 



5 



tivity ( |King et al. 2005[ l. The optical emission, on the other hand, 
is consistent with an external shock afterglow. 



4. 1 . X-ray flares 

The spectral evolution exhibited in the hardness ratio plot (fig- 
ureO supports internal shocks as the source of the flaring X-ray 
emission as suggested for GRB 050607 by Pagani et al. (2006). 
Though in this case the spectral index of the flaring period is 
similar to that of the late afterglow, the afterglow is softer, fol- 
lowing the trend exhibited by the hardness ratio plot. It should be 
noted also, that the spectral index of the flaring period is an aver- 
age over that period. The difference between the two values and 
the lack of spectral evolution in the late afterglow certainly sug- 
gests two distinct origins for energy emission. Combined with 
the overlap of the last y-ray and first X-ray flare, it seems likely 
that these two would share internal shocks as the common emis- 
sion mechanism. 

|Liang et al. (2006) show that X-ray flares ca n be modeled by 
the curvature effect (Kumar & Panaitescu 2000) which causes a 
decay of F oc (t - t & 0~^ 2+ ^ superimposed on the regular after- 
glow decay of F oc f~ a , where f e j is the energy injection time of 
that flare. We tested the flaring X-ray data against the curvature 
effect, using the average spectral index over that period and find 
ejection times for the 3 main X-ray flares at ~ 100 s, 200 s & 
380 s of 72 ± 6 s, 155+ 16s & 3 10 ± 16 s respectively. These are 



4.3. Early optical emission 



consistent with Liang et aL] - insofar as the energy ejection times 
are at the the start of the rising segment of the flare - supporting 
the claim that the flares are indeed due to internal shocks. 



4.2. Spectral indices & host extinction 

From the fitted flux of the FTN optical data points at 5000 s, 
corrected for galactic extinction and extinction due to Lyman 
absorption dMadau 1995K we find an optical spectral index of 
p opt = 3.1 + 0.4. Converting the X-ray flux at the same time 
to mJy using the X-ray spectral index fix - 1.14, gives a flux 
^i.732kev = 3.27 x 10~ 3 mJy. These fluxes correspond to an op- 
tical to X-ray spectral index, fi op t-x = 0.3 + 0.1, confirming the 
optical to X-ray flux ratio of F Vy0p JF Vt x ~ 10 of lStanek et al.l 

The large differences between fi opt , fix and fi opt -x imply that 
there could be an amount of host extinction that we have failed 
to take into account. To estimate this host extinction, we extrap- 
olate the X-ray flux to optical magnitudes to measure the optical 
extinction, A V (i+,) above that of Galactic. From this we can cal- 
culate the corresponding value of Eb v assuming SMC, LMC 
and Galactic extinction curves dPei 1992) in the host galaxy 
(Table[4). We do this in the limits that the cooling break is below 
the optical (fi pt-x = fix = 1.14) and that the cooling break is far 
above the optical but below the X-ray region (fi opt -x = fix~0.5 = 
0.64 below 10 17 Hz; |Zhang & Meszaros 2004) . The R & I band 
values of Eb v are consistent with each other at the 2<x level in 
all cases, so one case cannot be favoured over the other. In the 
limit where the cooling break is just below the X-ray region, we 
have lower limits on the host extinction. 

Using the host column densities found for the various extinc- 
tion curves (Table [2) we calculate the ratio of column density to 
optical extinction (Nh/Eb-v, Table |4) in the host and compare 
to the values of the SMC, LMC and Galaxy (IPei 1992) . The ex- 



|Stanek et al. (2007) suggest that the optical peak at ~ 540 s is the 
onset of the external shock - i.e. the deceleration time, when the 
jet has swept up enough of the circumburst medium so that the 
afterglow dominates emission - however if this was the case we 
would expect the X-ray to exhibit a peak at the same time. While 
the X-ray light curve at that time is highly obscured by flares, it 
does seem that the flares are superimposed on the already decay- 
ing X-ray afterglow. 

Since early time data has become more common due to the 
rapid dissemination of burst information from Swift, there has 
been much discussion regarding the reference time, fo, for the 
onset of the power-law decay of the afterglow. Most commonly 
it is taken as the trigger time To of the instrument which detected 
the burst, which is clearly instrument dependent. We expect that 
the light curve will evolve as F oc (t - toY", and that it will 
exhibit a break if it has been assigned the incorrect reference 
time ( |Quimby et al. 2006) . 

To test whether this is a plausible scenario, we fit the op- 
tical light curve at times > 540 s letting to as a free parame- 
ter. The best fit of this suggests that reference time is in fact 
the trigger time of the burst but since the x 2 function is very 
flat in that region (T <k 500 s), a definite minimum is difficult 
to estimate. This is in agreement with the theoretical prediction 
of Kobayashi & Zhang (2007 ) who claim that the reference time 
should be taken as Tq. The X-ray slope, which is at later times, 
is unaffected by changes in the reference time much less than the 
start time of the observations. 

If the deceleration time is, as we claim, close to 
the trigger time, the flat optical light curve from 70s - 
310 s is then likely due to external shock. The blast wave 



model (Zhang & Meszaros 2004 1 in a slow-cooling wind-driven 



pected ratios are 4.5, 2.4 & 0.48 x 10 cm respectively, clearly 
well below the measured values. This is in agreement with pre- 
vious work (Galama & Wijers 2001 ) suggesting dust destruction 
in the circumburst environment. 



medium predicts a temporal slope a opt = 0, in agreement with 
the observations (a opt = 0.10 + 0.10), for all values of p if 
v < v m < v c where v m is the peak frequency. 

We suggest that the rebrightening of a opt = -0.9 in the op- 
tical light curve after 310 seconds could be explained by the 
change from a wind-driven medium to one of a higher, constant 
density as detailed by |Pe'er & Wijers (2006) . If, as expected, the 
progenitor of this burst is a massive star, the circumburst medium 
is composed of the wind from the star moving into the ISM. This 
causes a forward shock into the ISM and a reverse shock into 
the wind, separated by a region of shocked wind with a constant 
density profile. When the blast wave crosses the reverse shock 
- shocked wind discontinuity, there is a drop in the cooling fre- 
quency. This may shift the optical out of the v < v m < v c regime 
into either v < v c < v m or v c < v < v m . If this is matched by the 
standard evolution of the peak frequency, the optical may end up 
in the v m < v < v c or v mc < v regime. The shocked wind has a 
higher density than the unshocked wind, necessitating a flux in- 
crease, causing the observed rebrightening. However, as we shall 
show in 94.51 the fine tuning of the energetics and circumburst 
medium parameters make this explanation of the rebrightening 
less likely. 



4.4. Late afterglow 

The X-ray spectral index of the late afterglow,/? = 1.14 + 0.03 
( 93.2b , implies an electron power-law index, p = 2.28 + 0.06, 
assuming that the cooling and peak frequencies are below the X- 
rays (Zhang & Meszaro s 2004) . If the observed X-ray break was 
a jet break, we would expect the temporal index after the break 
to be ax - p - 2.28 which is inconsistent with the observed 
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Table 4. Host extinction and column density to extinction ratio assuming SMC, LMC & Galactic extinction curves in the host as 
discussed in 34.21 In the case where v < v c , v c = 10 17 Hz. 
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decay of a — 1.35. This, in conjunction with the early time of 
the break (1.5^ - | hr in the rest frame of the GRB), makes a jet 
break highly unlikely. In the absence of a jet break, the electron 
power-law index, p, implies an X-ray temporal decay index of 
ax = 1.21 +0.05 in all media, which is not consistent with either 
of the broken power-law indices (Table 01. It could however, be 
in agreement with the single power-law. 

Table 5. Observed late afterglow temporal decays. The possible 
X-ray break is at 7.4*j - i hr. 



X-ray (power-law) 1.11 ±0.02 
X-ray (pre-break) 0.89 ± 0.05 
X-ray (post-break) 1.35 ± 0.06 
Optical 1.154 ± 0.013 



Assuming that the optical is in the same frequency regime 
as the X-rays i.e. above the cooling and peak frequencies, the 
expected optical temporal decay and spectral indices are pre- 
dicted to be the same as the X-ray indices, ax = a opt = 1.21, 
fix = ySopt = 1.14. This temporal decay index is in agreement 
with the observed value of a op t = 1.15 + 0.04 and consistent 
with the spectral index assuming host extinction. If the optical 
is not in the high-frequency regime, i.e. above either the cool- 
ing or peak frequencies, we can rule out a fast cooling regime 
(v c < v < v m ) as that would imply a = 0.25, which is clearly an 
underestimate. A slow cooling regime (v m < v < v c ) implies a 
decay of a — 0.96 + 0.05 in a homogeneous circumburst medium 
or one of a = 1.46 + 0.05 in a wind driven medium, each with a 
spectral index of ft = 0.64 ± 0.06. 

4.5. Blast wave energetics & circumburst medium 

From the measured quantaties we can constrain various physical 
parameters of the relativistic blast wave and its surroundings; the 
jet opening angle, the energy of the blast wave, the density and 
structure of the circumburst medium, the energy in the magnetic 
field and the relativistic electrons that emit the synchrotron ra- 
diation. In order to obtain these constraints we use the formulas 
from Van der Horst & Wijers (2007), while adopting p = 2.28 
and a luminosity distance c/l = 7.97 x 10 28 cm (£2 M = 0.27, 
Q A = 0.73, H = 71 kms- 1 Mpc" 1 ). 

Since the temporal behaviour of v m , v c and the peak flux 
F v ,mm depends on the structure of the circumburst medium, we 
derive constraints on the physical parameters in two cases: a 
homogeneous medium, in which the density is constant, and a 
stellar wind, in which the density as a function of radius is a 
power-law with index k = -2. In the homogeneous case F y>max is 
constant, and v m and v c are decaying power-laws in time, which 



gives the constraints: F v>max > F n , v m < v\ and v c < v\, all at 600 
seconds. A fourth constraint comes from the measured I-band 
flux at 600 seconds, corrected for galactic and host extinction, 
F Vl - 150 mJy. In the stellar wind case three of the four con- 
straints are the same, except for v c , because its value increases 
in time; this means that v c < v\ at 7 x 10 3 seconds. 

The limits on the isotropic equivalent energy we derive are 
lower limits of 5 x 10 54 erg for the homogeneous medium and 
2 x 10 55 erg for the stellar wind. Since these limits are quite 
large, we assume that the energies are not much larger. This re- 
sults in values for the fractional energies in electrons and mag- 
netic field of s e ~ eb ~ 10~ 2 (homogeneous) and e e ~ «b ~ 
5 x lO -3 (wind). The density of the homogeneous medium is 
then ~ 10 cm- 3 , while the wind density parameter A t ~ 1, i.e. a 
mass-loss rate of lO -5 solar masses per year and a wind veloc- 
ity of 10 3 km s . With these numbers and the lower limit on the 
jet-break time of ~ 10 6 s we derive a lower limit on the jet open- 
ing angle of 8° (homogeneous) or 2°(wind), resulting in collima- 
tion corrected energies of 5 x 10 52 and 1 x 10 52 erg respectively. 
These energies are quite large, but not impossible, especially if 
you consider that the early optical emission could be due to sus- 
tained energy injection by the progenitor, providing even more 
energy for the blast wave than there was for the prompt emis- 
sion. Furthermore, we note that a blast wave at this redshift and 
suffering significant host extinction has to be very energetic to 
produce the observed fluxes. 

We have suggested that the early optical behaviour could be 
explained by the transition from a massive stellar wind to a ho- 
mogeneous medium. To get the correct temporal slopes this im- 
plies that v\ < v m < v c from 70 - 3 10 s, and v m c < v\ after 540 s. 
This requires some fine-tuning of the parameters, in particular v c 
has to decrease significantly because of the density jump, and at 
the same time v m has to pass through the I band, since the peak 
frequency does not depend on the density and only mildly on the 
density structure. The isotropic equivalent energy that one de- 
rives in this case is even higher than before, namely 6 x 10 55 erg, 
while e e ~ 4 x 10~ 3 and sb ~ 3 x lO -3 . The density of the homo- 
geneous medium is again ~ 10 cm- 3 , but A„ ~ 0.15 in this case. 
From these parameters we get a lower limit for the jet opening 
angle of 6°and a collimation corrected energy of 3 x 10 53 erg. 
This value for the energy is very high and difficult to accom- 
modate within current progenitor models, making the transition 
from a wind to a homogeneous medium a less likely explanation 
for the early optical emission. 

We note that with these physical parameters the synchrotron 
self-absorption frequency, v a , lies well above 8.64 GHz at 3.95 
days when the 3cr limit was obtained by Frail (2006). The radio 
limit is hence consistent with both a homogeneous and wind- 
driven circumburst medium. 
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4.6. Energy injection 

Although the X-ray light curve is best fit with a broken power- 
law, it is difficult to accommodate the temporal slopes within 
our current afterglow models. The pre-break slope is shallower 
than the simultaneous optical slope, which is only possible if the 
circumburst medium is a stellar wind and v m < v opt < v c < vx- 
This would imply that the optical slope should be a op t = 1 .46 
and that the X-ray slope ax = 1.21, which is much higher than 
the observed pre-break values. Continued energy injection from 
the central engine (Nou sek et al. 2 006 ) would cause this slope 
to be shallower for the period of injection before reverting to 
the original value. For the optical temporal slope this means that 
the energy has to increase as f°- 37±0 - 07 , while the pre-break X-ray 
slope indicates /°- 30±008 . The break in the X-ray light curve at 
~ 7.4 hr is then the end of the energy injection phase. The post- 
break X-ray temporal slope is consistent with the spectral slope 
at a 2cr level for v m c < vx- The temporal break can not be due 
to the passage of v c , since there is no spectral change observed 
between the pre-break and post-break spectra. 

From the constraints on the ordering of the peak, cooling and 
observing frequencies we can again obtain limits on the physi- 
cal parameters of the blast wave and its surroundings. For the 
calculations we adopt that the isotropic equivalent energy is pro- 
portional to f° 34 , changing the temporal scaling laws of the peak 
flux, peak frequency and cooling frequency into F Vjmax oc f -0 ' 33 , 
v m oc f 67 and v c oc f -1,33 . The values for the physical parame- 
ters we derive are quite similar to the values in 34.51 sb ~ 10" 3 , 
e e ~ 10~ 2 , A* ~ 1, the isotropic equivalent energy at the end of 
the energy injection phase ~ 5 x 10 55 erg, the jet-opening angle 
> 2°, and the collimation corrected energy ~ 2 x 10 52 erg. So 
again, while the energy requirements are at the high end of the 
distribution, they are not unprecedented. 

5. Conclusion 

We have analysed the optical, X-ray and y-ray data of 
GRB 060210 from the time of the prompt emission up until the 
X-ray afterglow was no longer detectable by Swift-XRT. 

The early flaring X-ray & y-ray emission implies internal 
shocks, indicative of prolonged engine activity. The light curves 
of the two regimes match well, even displaying a common flare 
at ~ 100 s. The agreement of the X-ray flares to the curvature 
effect and the clearly demonstrated spectral evolution support 
the claim that the flares are indeed due to internal shocks. The 
simultaneous optical emission, on the other hand, is consistent 
with an external shock afterglow with a frequency below that of 
both the peak & cooling frequencies. We have shown that such 
an external shock's deceleration time is close to the trigger time 
and should hence dominate the optical emission at the time of 
our observations. 

The late afterglow (> 600 s) may be explained by either of 
two models with an electron power-law index, p = 2.28, derived 
from the X-ray spectral index, fix'- a single temporal power-law 
in both optical and X-rays with both observing bands above the 
cooling frequency; or a broken power-law in X-rays in which 
optical and pre-break X-ray temporal slopes can be explained 
by continued activity of the central engine. In the latter case the 
cooling frequency is between the optical and X-ray frequencies 
and the observed X-ray temporal break indicates the termination 
of the central engine activity. Though both are consistent with 
the temporal and spectral slopes, we favour the continued central 
engine activity since the x 1 of the single power-law fit is signifi- 
cantly worse. The collimation corrected energy requirements of 



~ 10 erg in both cases, are at the high end of the distribution 
for GRBs but are certainly not unprecedented. 

Comparing the column density in the host to the optical ex- 
tinction we find a higher than expected value which supports the 
notion that GRBs can cause destruction of dust in the circum- 
burst environment. 
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